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The paper examines part of the ground state diagram of the extended Hubbard model, with the 
on-site attraction U < and intersite repulsion W > in the presence of charge density waves, 
' superconducting and ?7-superconducting order parameters. We show the possibility of the stabiliza- 

, tion of the mixed state, with all three nonzero order parameters, in the model with nearest neighbor 

■ interactions. The other result of the paper is application of the exact solution of the Schrodinger 

04 ' equation for the two-electron bound state, as an additional bound for the phase diagram of the 

j^jr^, model, resulting in the partial suppression of the superconducting state of the s-wave symmetry, in 

. favor of the normal state phase. 
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PACS numbers: 71.10.Fd, 71.10.Hf, 71.27.+a 



INTRODUCTION 



The Hubbard model is known for its universality. The model and its extensions have been used in research of 
magnetism, superconductivity and especially high temperature superconductivity (HTS), other various phases in 
^ ■ the solid, charge orderings, phase separation etc., in materials like high temperature superconductors, bismuthates, 
Chevrel phases, amorphous semiconductors, heavy fermion materials, systems with alternating valence to name the 
few (for a review see, e.g., Ref.Q). Despite intensive research the results concerning the model come mostly from 
different approximations or numerical simulations. The exact results are scarce and concern mostly specific situation 
''O of one dimension (ID) or specific band fillings or specific limits: zero bandwidth, infinite U or infinite d. We can 
Ch mention the exact solution in ID obtained by Bethe ansatzQ, Nagaoka ferromagnetism in repulsive, half-filled systems 
^ , with one hole'^, flat band ferromagnetism Q, Lieb's ferrimagnetisml^, some bounds on correlation functions'^, 

and a statement of necessity of bound states for existence of extended s-wave superconductivity in 2D systems with 
low electron density n, by Randeria[^. Let's also note that the mean-field ECS equations for superconductivity in 
' ^ ' the Hubbard model with effective attractive interaction between electrons, in the limit of vanishing electron density 
, turn into Schrodinger equations, which can also be solved exactly ^3 ■ 

' An interesting, unresolved question concerning Hubbard model is: what happens, when we increase the strength of 
\^ the interaction? What phases do we obtain? Or does the system prefer phase separation? That was the conclusion 
OO of the paper by Robaszkiewicz and Pawlowski llj. By the broken symmetry Hartree-Fock approach they established 
that the attractive extended Hubbard model (with nearest neighbor - nn - density-density interactions of strength 
W added) in two dimensions for arbitrary electron density and moderate values of \W\ is superconducting in the 
ground state but becomes phase separated with increasing \W\ (except from n = 1, where the ground state is in 
' charge density wave phase for any > 0) . In the case oiW<0 they obtained the phase of electron droplets and for 
^ . W > a, phase separation between charge density waves (CDW) with n = 1 and a singlet on-site superconductivity 
I ' (SS) (extended s-wave was not considered). The free energy of this phase separated state (PS[CDW/SS]) turned out 
^ , to be lower than the energy of the homogenous mixed state (M(CDW-I-SS)). Stabilization of M state required longer 
than nn repulsive intersite interactions (W2) or increasing \U\. For U = —00 PS and M states were degenerated even 
Q . for W2 = 0. 

^ ' In this paper it will be shown that including 77-pairing (superconducting pairing with singlet pairs with center-of- 
, mass momentum 7r)|l2L ll^l Il4l| into the Hamiltonian stabilizes mixed phase, already with intersite interactions W 
restricted to the nn (in disagreement with Rcf. llj). It turns out that the phase separation PS[CDW/SS] exists only 
^ , below certain critical value of W, above which the M phase can be stabilized. 



LOW DENSITY LIMIT 



It can be shown, that in the low density limit the BCS equations go over to the Schrodinger equation for the bound 
pair^^. In the center-of mass coordinate system the wave function of the two particle bound state can be expanded 
as the sum of states with the definite relative momenta. Using that expansion we can obtain a set of self-consistent 
equations for the pair wave function in a position space, in terms of lattice Green functions [H. Il5l| . For the hypercubic 
lattices these equations were solved (see, e.g., |l|) and it has been found out that in one and two dimensions pairs 
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for W = bind for any negative U, while in three dimensions there is critical value for W. In the cited reference the 
formula was given for the Wcr as a function of U : 



\Wcr 

2t 



2D 

IT 



+ {c-ir' (1) 



where C — ^/Nj^ki^ ~ lk/z)~^ is the Watson integral. In Ref.0] this formula was given only for the case of < 
but in fact it is valid for any U and W . What's more, there exists analogic formula, applicable not only in the low 
density limit but for any electron density. It was derived from the Hamiltonian after Hartree-Fock linearization so it 
is not exact, but goes over into the exact Eq. for n 0. The formula is given by: 



Wcr^— , — TTT^TTT' where ^ = / xp{x) dx (2) 



Ai(l-n) + 8</-2^VC/' 
or for the rectangular DOSlij: Wcr{l + {n - 1)2(1 + 16t/U)) = -it. 

PHASE SEPARATION 

The energy of a phase separated state, Eps, is calculated according to the formula: 

Eps ^mE+{n+) + {l-m)E_{n_) (3) 

where E±{n±) are the values of < H > /N of the two separating phases at the n = n± corresponding to the lowest 
energy homogenous solution for a given phase, m and (1 — to) are fractions of the system occupied by the phases E^ 
and E^ respectively. There is also a similar equation connecting charge densities in the two coexisting phases: 

TOn_|_ + (1 — m)n_ = n (4) 

A very useful is also Maxwell construction: 

9£;_(n_) E+{n+) - E^{n^) 



(5) 



which is equivalent to the condition of minimization of Eps respectively n- with fixed 7i+ and 7i. Taking into account 
Eqs l@J and ©, Eq. Q can be rewritten as: 

Eps = E+{n+) - (n+ - n)^^^^ (6) 

There is also another, equivalent but more intuitive approach using Gibbs energies G = F — fin. The independent 
variable for G is chemical potential /i instead of charge density n for F. Let's note, that instead of the analogue of 
Eq. I0J, with fi± instead of n±, we have /i+ = /i_, because this is the condition of equilibrium, or rather the definition 
of equilibrium! (Eq. does not apply also because /i is not an extensive quantity). There is only one value of 
chemical potential, fixed and the same for both phases. This way the whole construction simplifies. The procedure 
is as follows: we find solutions for the homogenous phases, then we plot their Gibbs energies vs fj, in one figure and 
we look for the intersections of G"s belonging to the different phases. Such intersection means a coexistence of the 
two phases with the same /i, i.e., in equilibrium. The same value of /i usually translates into different values of n+ 
and 7i_ in the two phases under consideration and that's nothing else but the (definition of) phase separation, which 
appears on the F-n phase diagram for the electron density n, constructed according to the Eq. Q). 

Some remarks are due: 

— let's note that the use of G and chemical potential makes it clear that 7i+ and 7i„ are independent of n, as there is 
usually one intersection point of Gibbs energies. In other words, tti is a linear function of n: 

TO = ri (7) 

n_|_ — n_ — 7i_ 
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— let's note that Eq. 10) gives boundaries for Eps- By the definition in Eq. (jSJ, Eps must be somewhere between 

£'+(n+) and E-{n-). By the same reasoning that led us to the Eq. ^ we can say that the dependence of phase 
separated states on n should be strictly linear, with the same coefficients as in the Eq. {71). That's nothing else 
but the Maxwell condition restated 

— the same considerations as above apply to the electron density: phase separated state can exist only for charge 

densities n_ < n < This is an important bound, which can be easily overlooked, as there are formal 

solutions for free energy of phase separated state also outside that range 

— a simple implication of the preceding points, which we will call Theorem 1: we can not obtain a phase transition by 

changing n between two phase separated states having "a state in common" for the same fixed density n+, i.e. 
we can not have phase transition between states like PSl[CDW(n=l)/phasel] and PS2[CDW(n=l)/phase2]. 
The lines depicting Fpsi and Fps2 will begin in the same point and will extend in different directions. 

— and last but not the least is the notion, that we have also bounds for homogenous phases: the phases can exist only 

for the positive compressibility, i.e. in the range of /i, for which dfj,/dn > 

THE MODEL 

We start from the extended Hubbard model: 

ij(T i ij i 

in standard notation. We use Hartree-Fock approach. In the mean-field Hamiltonian we retain charge density order 
parameter Ac and singlet superconductivity order parameters: Aq and Ag, for on-site and 77-pairing respectively (see 
also Ref.[l3): 



Ho = "^{Sk- ^l)cl^Cka + ^^i^ccl„Ck+Qa+h.C.) + 

k(7 ka 

^(Aocl-fcl,^ + h.c) + ^(Agci+Q-fcl,^ + h.c) + C (9) 



k 

where the order parameters are given by: 

A, 



l^(-zW^+^)<4+Q^Cfe. > (10) 

k(7 

^0 = -^Xl < ^-''^^iCfe^T > (11) 



k2 



= -^"^ < ^k2lC-k2+Q1 > (12) 
fe2 

where z is number of nearest neighbors, JI — ji — + zW)n, Wk = W^7fc, £k — —tlk, Ik = '^Yl'L cos/cqi 01 G [x^y, z) 
and d ~ z/2 is dimensionality of the hypercubic lattice; £/j G< — eoi^o > where eg = D = zt. In the following all 
values of dimension energy will be given in units of D and the simplest, zero temperature forms of the equations will 
be most often used. The constant C takes the form: 

C=ku + 2WzW - _ ^ _ ^ (13) 

4^ ' U-2Wz U U ^ ' 

After diagonalization we obtain Hamiltonian Hq with two eigenenergies A^: 



At = V ^fe + + + A| + A2 ± 2^el{jl^ + A|) + (pAe - AoAq)^ (14) 
In terms of the eigenstates of Hq the free energy is given by the standard formula: 

Fo = lnTrexp(-/3Ho) + /iiVe (15) 
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what yields: 

^=-pin-l) + C-^Y. 1^(4 cosh ^ cosh ^) (16) 

The equations for the order parameters and chemical potential are obtained in a standard way, by minimization of the 
above expression. For the most general mixed, homogenous phase with all the order parameters present the equations 
take the form: 

2{n - 1) = -lilo + J2) + Ac(mAc - AoAq) Jo (17) 
"4Aq 



U 
-4Ao 

U 
-4A, 



U -2Wz 



AQ{Io + J2)-AoiTiA,-AoAQ)Jo (18) 
Ao/o-Aq(7IA, -AoAq)Jo (19) 
Ae/o+7i(MAc- AoAq)Jo (20) 



where 



/n = 



-y(- L) ^fc (22) 

^4(m^ + A^) + (7IA,-AoAq)2 



where are given by Eq. (|14|l . The equations for the mixed phases with smaller number of constituent phases (i.e., 
M(CDW+SS), M(CDW+SQ), M(SS+SQ)) as well as for pure phases can be obtained from Eqs. by just 

dropping out respective order parameters. For the detailed form of the pure phase equations the reader is referred, 
e.g. to Ref. 
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RESULTS 



The results were obtained using rectangular density of states (DOS) and z = 4 as an approximation for two 
dimensional square lattice. In Figs 1 (a) and (b) we plotted the results for the thermodynamical potentials. In (a) 
we have shown all the solutions for Gibbs energy (vs chemical potential), without applying yet the bonds for n± or 
checking the sign of dfi/dn. All homogenous phases are shown, except from Ty-pairing (denoted by SQ) and mixed 
M(SS+SQ) (SS denoting singlet, on-site superconductivity), which do not have solutions for U > —2. Nevertheless 
the phases: M(CDW+SQ), M(CDW+SS), upper branch of M(CDW+SS+SQ) (see lower inset) and upper branch of 
CDW for 71 ^ 1 must be excluded, because of the negative compressibility, connected with the downward curvature of 
these lines (compare also the upper inset - decreasing and increasing part of the plot corresponds to the opposite signs 
of compressibility and to the two branches of M(CDW-|-SS+SQ), from now on denoted as M3, in the main figure). 

The remaining curves (and CDW one, for comparison) are a basis for plotting the free energy F-n diagram, shown 
in (b). The plot contains curves in the proper range of n, with all the bounds already applied. Apart from all the 
homogenous phases fulfilling the condition dpi/dn > (what corresponds to the downward curvature for F{n)) we 
can see also phase separated states indicated in (a) by Gibbs energy plots crossings: PS[CDW/SS], PS[CDW/NO], 
PS[M3/SS] (and PS[NO/NO], of slightly different nature). The states phase separated with CDW(n=l) are connecting 
the point of energy E'~^^^ {n = 1) with respective curves of the other state coming into separation, in a tangential 
way, in accordance with Maxwell construction and our earlier remarks. Let's note, that the CDW inside PS states is a 
phase with n = 1 while CDW in mixed states is for n 7^ 1. The states of the lowest energy are: SS, M3, PS[CDW/SS] 
and PS[M3/SS] and it will be so also for other values of parameters. In the inset there are the energy differences 
between the states: PS[CDW/SS], PS[M3/SS] and M3, difficult to discern in the main figure. We can see that the 
phase separated state PS[M3/SS] has the lowest energy of all the three (if exists). Its appearance is an unavoidable 
consequence of Theorem 1 in the case when the phase M3 starts to win with PS[CDW/SS]. The phase diagrams are 
shown in Fig. 2 (a) and (b). 

For not too large \U\ and W the ground state is in PS[CDW/SS] or SS state, depending on n, in accordance with 
Ref.^3|. When we increase W we reach a threshold, smallest for the half-filled band, above which the M3 phase 



5 




FIG. 1: (a) Gibbs energies G vs bare chemical potential ^, for U = —1 and W = 0.2. Lines are described in the figure. In 
the lower, left-hand corner in the inset, there is an enlargement of the area where the curves cross. In the inset in the upper 
right-hand corner there is plotted the dependence of chemical potential ^ on charge density n in M3 phase, (b) The free energy 
F vs electron density n for U = —1 and W — 0.2. Lines are described in the figure. In the inset there are the differences of 
free energies plotted - full line: F{PS[CDW/SS]) - F{M3), dashed line: F{PS[CDW/SS]) - F(PS\MZl SS\), dotted line: 
F{M3) — F{PS[M3/SS]) for the range of n, where dfi/dn > for each phase. Horizontal dotted line is the just a guide for an 
eye. 



is stabilized. Appearing of the M3 phase is connected with replacing PS[CDW/SS] by PS[M3/SS] phase, what is 
required by the Maxwell construction (see Eq. (|SJ) and Theorem 1). With increasing \U\ the threshold and the extent 
of PS[CDW/SS] phase are decreased and, for large enough \U\, the M3 phase is stabihzed for all > at n = 1, as 
shown in Fig. (2b). The area of PS[CDW/SS] disappears at all, in favor of mixed phase. 

In both Figs (2a) and (2b) there is also a line "cutting off' the upper right-hand corner. This is the boundary 
for existence of bound pairs of extended s-wave symmetry, given by Eq. In view of Randeria's notion^, about 
necessity of such pairs for existence of superconductivity of this symmetry it can be used as an additional bound for 
the phase diagram Above these lines no superconductivity can survive so the NO phase enters diagram (to be 




FIG. 2: 

(b) C/ = 
phases. 



Ground state phase diagram of model Eq. ||UJ calculated for the rectangular density of states for (a) U — —0.925 and 
-L2. For the n = 1 the ground state is CDW phase, a smooth limit of both M(CDW+SS-hSQ) and PS[CDW/SS] 
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sure of that we must consult again the free energy diagrams). Let's note, that it is only an approximation, as we 
did not calculated extended s-wave superconductivity but only on-site, pure s. Preliminary results show nevertheless 
that the parameter of extended superconductivity is two orders of magnitude smaller than Aq along the boundary 
PS/SS state. Also 77-pairing can have extended component, which will be stabilized by 14^ > 0, due to negative sign 
at the extended part of the pairing potential, unlike the case of extended s-wave superconductivity with Q = 0. This 
fact will still enhance the main tendency outlined in this paper. 

Inclusion of the normal Fock parameter into the calculations will bring about the increase of the bandwidth, thus 
effectively pushing the system into weak-coupling limit. It can cause the increasing of thresholds but should not 
change qualitatively the main results of this paper. 

CONCLUSIONS 

In the paper a phase diagram of the extended Hubbard model with attractive on-site and repulsive, nn, intersite 
interactions was analyzed in the Hartree-Fock approximation. A possibility of stabilizing the mixed phase in the 
model with short range, nn interactions only, was shown. Stabilization requires the existence of the three orderings: 
superconductivity, 77-superconductivity and charge density waves. The extent of the mixed phase increases with 
increasing \U\. 

The other result of the paper comes from the solution of the Schrodinger equation for bound pairs of extended 
s-wave symmetry . The area of the normal phase is introduced into the phase diagram due to this criterion. 
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